Compelling evidence of renormalons in QCD from high order perturbative expansions 
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We compute the static self-energy of SU(3) gauge theory in four spacetime dimensions to order 
in the strong coupling constant a. We employ lattice regularization to enable a numerical 
simulation within the framework of stochastic perturbation theory. We find perfect agreement with 
the factorial growth of high order coefficients predicted by the conjectured renormalon picture based 
on the operator product expansion. 
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Little is known about properties of quantum field the- 
ories from first principles. This is particularly so for 
asymptotically free gauge theories such as quantum gluo- 
dynamics. One of the most salient features of this theory 
is the confinement of charged objects. Yet this prop- 
erty has not been proven, and the best evidence comes 
from the linearly rising static potential at large distances 
obtained in lattice simulations. Another expected prop- 
erty is the asymptotic nature of perturbative weak cou- 
pling expansions. In four dimensional non-Abelian gauge 
theories one particular pattern of asymptotic divergence 
should be determined by the structure of the operator 
product expansion (OPE). It is usually named renor- 
malon m or, more specifically, infrared renormalon. Its 
existence has also not been proven but only tested as- 
suming the dominance of /3o-terms, which amounts to 
an effective Abelianization of the theory, or in the two 
dimensional O(A^) model [H, where it is suppressed by 
powers of 1/N. Moreover, the possible non-existence or 
irrelevance of renormalons in Quantum Chromodynamics 
has been suggested in several papers, see, e.g. [l, 0] and 
references therein. This has motivated dedicated high 
order perturbative expansions of the plaquette, e.g. [5|- 
@], in lattice regularization, with conflicting conclusions. 
Powers as high as a^'^ were achieved in the most recent 
simulation |9j. However, the expected asymptotic be- 
haviour was not seen. If confirmed, this non-observation 
would cast doubt on the well-accepted lore of the OPE 
and renormalon physics (see [lo| for a comprehensive re- 
view), and would significantly affect the phenomenologi- 
cal analysis of data from high energy physics experiments 
on the decay of heavy hadrons, heavy quark masses, the 
running coupling parameter, parton distributions, etc.. 
Therefore, this issue should be clarified unambiguously. 

In this letter we present compelling numerical evidence 
that the expected renormalons indeed exist not only in 
models but in real gluodynamics. We also argue why 
previous analyses based on the plaquette have failed to 
detect them. The vital and new ingredients of our study 
are as follows. 

(a) We consider a perturbative series whose leading 
renormalon is dictated by a dimension d — 1 operator. 



rather than by the d = 4 plaquette. 

(b) Using a higher order integrator and employing 
twisted boundary conditions, among other improve- 
ments, we are able to obtain results of unprecedented 
precision on an extensive set of spacetime volumes. 

(c) We carefully extrapolate to the infinite volume limit, 
thoroughly investigating finite size effects. 

Perturbative expansions in powers of a, 
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are believed to be asymptotic and not Borel summable 
in QCD, due to the existence of singularities in the Borel 
plane (renormalons). Typically fc„ will diverge like aJ^nl^ 
with a constant ad- This divergence pattern of fc„ should 
not be arbitrary but consistent with the OPE associ- 
ated to a physical observable. Even though this factorial 
growth was originally discovered analyzing the Feynman 
diagrams that contribute to the large /3o approximation, 
the correct divergent structure can only be inferred by 
assuming that the perturbative series is asymptotic and 
complies with the OPE. The OPE fixes the positions and 
the structure of the renormalon singularities in the com- 
plex Borel plane, resulting in a more intricate pattern 
that cannot be obtained from the large /3o approxima- 
tion alone. Successive contributions fc„Q;" decrease for 
small orders n down to a minimum at no ~ l/darfja). 
Higher order contributions should be neglected and in- 
troduce an ambiguity of the order of this minimum term, 
fcrio""" ^ exp[-l/(|arf|Q;)]. 

Within the OPE an observable R{q, A) can be factor- 
ized into short distance Wilson coefficients Ci{q, i-i) and 
non-perturbative matrix elements {Oi{fj,,A)) of dimen- 
sion i: 



R = Co{q, fi){Oo{fi, A)) + Cd{q, /i)(Od(M, A)) 

\ q / 

(2) 

q, A and fi denote a perturbative, low momentum and 
factorization scale, respectively, so that g ^ /i ^ A. 
For the plaquette, (Oq) — 1 and the next higher non- 
vanishing operator is the dimension d — 4 gluon con- 
densate. In this case, the perturbative expansion of Co 
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cannot be more accurate than 0{A!^ /q^) which is exactly 



of the size of the fc„„ term since 
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with /3o — 11. The so-cahed leading infrared renormalon 
of this expansion cancels the ultraviolet ambiguity of the 
next order non-perturbative matrix element so that the 
physical observable R is well-defined. 

From this discussion it is evident that we should study 
series expansions with the smallest possible tiq or, equiv- 
alently, d. For d = 1 the perturbative expansion should 
start to diverge at an order no that amounts to about 
one fourth of that for the plaquette. This applies to the 
pole mass (see 11, l^]) and to the associated self-energy 
of a static source, which we consider here. The latter 
does not have a continuum limit, as it linearly depends 
on the ultraviolet regulator. Here we consider lattice reg- 
ularization with the Wilson gauge action [l3| and write 
the self-energy in the following way: 
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n>0 



a~^, the inverse lattice spacing, provides the ultravio- 
let cut-off. The large n behaviour of the coefficients c„ is 
regulator independent, universal and equal to the asymp- 
totic behaviour of the pole mass up to ©(e"^/") terms 
(due to subleading renormalons): 
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[n + b) 
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The coefficients b and si were computed in 14]. They 
read (see [l5| for details) 
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For a static source in the fundamental (triplet) repre- 
sentation the normalization constant N.^ is exactly the 
same as for the leading renormalon of a heavy quark pole 
mass. This renormalon is also related to a renormalon 
of the singlet static potential since these contributions 
cancel from the energy -E(r) = 2™ -f V{r) [l6l - [l8j |. For 
adjoint sources it corresponds to a specific combination of 
pole mass and adjoint static potential renormalons [l9| . 
The factor Nm is cancelled in the ratios 
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We obtain the expansion coefficients c„ of the static 
energy from the temporal Polyakov line on hypercubic 



TABLE I. Lattice geometries. Volumes with boldface time 
extents are expanded up to ©(a^"), the others up to 0{a^^). 
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lattices. We investigate volumes of Nt lattice points in 
the time direction and spatial extents of Ns points. For- 
mally we may introduce an anisotropy at ^ as- lu this 
case the lattice action, that is invariant under time or 
parity reversal, agrees with the continuum action up to 
0{a\,a\^ terms. The temporal and spatial lattice ex- 
tents in physical units are given by atNx and agNs, re- 
spectively, so that the only dimensionless combinations 
consistent with the leading order lattice artefacts are 
a1/{atNTf = l/iV| and l/7Vj. Therefore, within per- 
turbation theory, where we cannot dynamically generate 
additional scales, the leading order lattice artefacts are 
indistinguishable from 0{1/N^, 1/Ng) finite size effects. 

We choose periodic boundary conditions in time and, 
to eliminate zero modes and to improve the numerical 
stability, twisted boundary conditions [20l - |23j | in all spa- 
tial directions. The Polyakov line is defined by 
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where U^{n) « e'aKl('^+'^/^)''] e SU(3) denotes a gauge 
link in representation _R, connecting the sites n and n-\-ji, 
TLi e {0, . . . ,iVs-l}, 714 G {0, . . . , A'^r-l} and g = VAttcx. 
We implement triplet and octet representations R of di- 
mensions dn = 3 and 8. The link U4^(n) appears within 
the covariant derivative of the static action ■0D4'0, the 
discretization of which is not unique. We use singly 
stout-smeared (smearing parameter p — 1/6) covari- 
ant transporters instead of Ui(n) as a second, alternative 
choice, to demonstrate the universality of our findings. 

We remark that neither the lattice spacing nor the 
strong coupling parameter a enter our simulations 
explicitly. Numerical stochastic perturbation theory 
(NSPT) [25l-[27l| enables us to directly calculate coeffi- 
cients of perturbative expansions. We employ the vari- 
ant of the Langevin algorithm introduced in [2^ that 
only quadratically depends on a time step At. Extrap- 
olations to At = were performed on a subset of lattice 
volumes where we found agreement within statistical er- 
rors between all our extrapolated expansion coefficients 
and those obtained at Ar = 0.05. For the geometries 
listed in Table U we restrict ourselves to this fixed value, 
which, within errors, effectively corresponds to At = 0. 
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We expand the logarithm of the smeared and un- 
smeared Polyakov lines in different representations to ob- 
tain the corresponding static energies: 



P{Ns,Nt) = - 



\n{L{Ns,NT)) Ns.Nt- 
aNr 



(9) 



Fortunately, the dependence of this logarithm on Nt and 
Ns can be deduced and only a few parameters need to 
be fitted at each order: 
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^ [c„ + AW(iVs) + Ai2)(iVs,iVT)]a"+V') - 

n>0 



Ai^)=^^'n-^te"'+logs(:)(^5) 
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The logs^f''(A^5) are polynomials of ln(A^5) of order n — 1, 
with coefficient functions that depend on /-,- and the /3- 
function coefficients /3j where j < n—l. These terms are 
entirely determined by the renormalization group run- 
ning of a. The logs^^^"'^(A^5) are obtained in the same 
way. In the Nt — >■ oo limit An"^ is the dominant cor- 
rection while A^^' includes the leading 0{1/N!^,l/Ng) 
lattice artifacts discussed above. 

The term A^^^ originates from interactions with mirror 
images, see also [29|. This effectively produces a static 
potential between charges separated at distances aNs, 
but without self- energies. Therefore, we expect the high 
order behaviour of /„ and c„ to be dominated by one 
and the same renormalon. This can also be illustrated 
considering the leading dressed gluon propagator D{k) cx 
where k4 = 0. With the (formal) ultraviolet cut-off 
1/a and an infrared cut-off l/{aNs) this can be written 
as (ignoring lattice corrections). 



Pcx / dkk'^D{k) 

Jl/(aNs) 



(11) 



after perturbatively expanding D{k). When re- 
expressing a{{aNs)^^) in terms of a{a^^) we may con- 
sider two situations: 

(a) Ns > e". In this limit the last term of Eq. (fTTj) 
is exponentially suppressed in n and the renormalon can 
directly be obtained from a large order expansion of aP. 



(b) Ns < e" . The last term of Eq. (ITT|) is important and 
the renormalon cancels order-by-order in n. 

In present-day numerical simulations Ns < e", and the 
term A^^-* needs to be taken into account, in combination 
with c„. A similar phenomenon was numerically observed 
for the static singlet energy E{r) — 2m + V{r) 3, 30[. 
This teaches us that to correctly identify the renormalon 
structure of 5m, it is compulsory to incorporate the \/Ns 
corrections. So far, in studies of high order perturbative 
expansions of the plaquette the corresponding finite size 
terms have been neglected. Our fits indeed yield /„ ~ c„ 
for large n, in clear support of the renormalon dominance 
picture. 

In the lattice scheme /3o, Pi and /32 are known [sij. 
The effects of higher jSj start at O(a^), but this uncer- 
tainty in our parametrization quickly becomes negligible 
at high orders where the coefficients fj , governed by the 
d — 1 renormalon, will dominate. This can be quanti- 
fied systematically in a large n analysis [s^ . where any 
possible renormalon of the lattice /3-function is sublead- 
ing (d > 1). To check this assumption and to justify the 
truncation at (S^ we have performed fits including j3j for 
j < 0, 1 and 2 (see below). 

Starting at C(a^), one may expect additional finite size 
terms oc hi{NT / N s) / N s from a possible mixing of the an- 
titriplet interaction between mirror charges with sextet 
and higher representations, mediated by ultrasoft gluons, 
in analogy to the mixing of singlet and octet static poten- 
tials in potential nonrelativistic QCD (pNRQCD) f33| . 
These terms are subleading from the renormalon point 
of view [d = 3). Moreover, aNs provides an infrared 
cut-off to gluon momenta so that one would only expect 
such contributions in the limit Ns ^ Nt that we do not 
investigate and, indeed, we see no numerical evidence of 
them. 

Our data are sensitive to the 1/N^ correction terms 
within A^^^. However, including w„ or //i""' as additional 
fit parameters did not significantly improve the x'^-value 
and so we decided to omit the l/iV| and l/iV| terms. 
Note that these contributions, if present, can numerically 
easily be distinguished from 1 /Ns and become irrelevant 
at relatively small Ns, unlike logs(Afs)/^s terms. 

As a cross-check we calculate diagrammatically, 

Co 2.1172743570834807985970 . . . , (12) 
ci = 11.1425(25) , /o = 0.76962563284(2) , (13) 



/(^"^ = 0.14932(3) , wo - 
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0, 
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for the unsmeared Polyakov line. In this case /o™\ the 
l/iV| coefficient, does not vanish but it is small. For 
fundamental sources, Cq and ci were known diagrammat- 
ically before and C2 numerically 2^, 34 1. Our fit repro- 
duces these values. For adjoint sources the above coeffi- 
cients need to be multiplied by the factor Ca/Cf = 9/4. 

We exemplify the result of our global fit to the un- 
smeared triplet data obtained on all our geometries (see 
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FIG. 1. Comparison between the global fit and data for n = 9. 
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FIG. 2. The ratio c„/(nc„_i) for the smeared and unsmeared 
fundamental static self-energies, compared to the prediction 
Eq. (O at different orders of the 1/n expansion. 



Table m and orders of perturbation theory (with four pa- 
rameters per order) in Fig. [1] where a comparison to the 
n = 9 data is shown. We find smeared and unsmeared 
data to be well described by the fits, with reasonable 
X^/iVDF ~ 1-29 and 1.46, respectively. Note that the fac- 
torial growth found (and expected) for the coefficients /„ 
produces very sizable l/iVg terms at high orders. 

In Fig. [5] we compare the infinite volume extrapolated 
ratios c„/(nc„_i) to the theoretical prediction, Eq. ([T]). 
LO, NLO and NNLO refer to this prediction, truncated 
at 0(1), 0{l/n) and ©(l/n^), respectively. The data are 
robust to subtracting lattice artefacts (the A^^^ terms 
of Eq. ([in])) or to truncating at different orders in 
see Fig. [31 Particularly reassuring is the universality 
of the result; fits to smeared and unsmeared Polyakov 
loop expansions give the same large n behaviour, fully 
consistent with the dominance and universality of the 
infrared renormalon; smearing only affects the ultravi- 
olet behaviour. Fits to the octet representation data 
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FIG. 3. The same as Fig.[5]for the unsmeared data, truncating 
at different orders in f]j. 



also show exactly the same behaviour, again in agree- 
ment with the renormalon dominance picture. Also note 
that NSPT data for different orders are statistically cor- 
related. These correlations work in our favour. We post- 
pone the details of this to ^1] . 

Finally we determine the normalization of the pole 
mass renormalon, see Eq. ([SJ, and obtain N}^^ = 
18.6(4) for the smeared and N^f = 19.0(3) for the un- 
smeared static action. Converting this to the modi- 
fied minimal subtraction (MS) scheme, we find = 
AiatA^mV^MS ^ 0.65(2). This agrees remarkably well 
with the estimate ~ 0.62 of [13, Hlj from a MS 

scheme expansion up to 0{a^), in support of the claim 
that renormalon dominance starts at much lower orders 
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in the MS scheme (see for instance 
inary results from directly converting our lattice data to 
the MS scheme further reinforce this claim. Irrespective 
of the scheme, the heavy quark pole mass can only be de- 
fined up to an ambiguity of ~ 0.65 A^. A more detailed 
analysis is in preparation (s^ . 

In conclusion, we have obtained the static self-energy 
of SU(3) gauge theory in four spacetime dimensions to 
0{a^'^) in the lattice scheme. For orders rt > 9 we find 
perfect agreement with the factorial growth of the coef- 
ficients, as predicted by the conjectured renormalon pic- 
ture based on the operator product expansion. Note that 
this implies that, in the lattice scheme, we expect the 
renormalon dominance of the plaquette to set in at val- 
ues n > 4 X 9 that so far have not been realized in the 
literature. 
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